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Abstract. We demonstrate the occurrence of free quasi-particle excitations obeying 
generalised exclusion statistics in a BCS model with asymmetric pair scattering. The 
'^ ' results are derived from an exact solution of the Hamiltonian, which was obtained via 

^ , the algebraic Bethe ansatz utilising the representation theory of an underlying Yangian 

algebra. The free quasi-particle excitations are associated to highest-weight states of 
the Yangian algebra, corresponding to a class of analytic solutions of the Bethe ansatz 
equations. 
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'nI" . 1. Introduction 
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In a recent Fast Track Communication [T3] we discussed the phenonemon of generalised 
exclusion statistics [1011111112] in a non-Hermitian BCS model. The non-Hermitian 
K^ I Hamiltonian results from an asymmetry in the Cooper pair scattering amplitudes, 

depending on whether the scattering is from high energy to low energy states or vice 
versa. As a result of being non-Hermitian, the Hamiltonian admits complex eigenvalues 
for particular choices of the coupling parameters. However as the Hamiltonian is a 
real-valued operator non-real eigenvalues, when they occur, arise as complex conjugate 
pairs. In this sense the spectrum of the Hamiltonian mirrors the property of admitting 
an unbroken "PT-symmetric phase and a broken PT-symmetric phase which has been 
discussed in multiple other models, e.g. see articles contained within the Journal of 
Physics A Special Issue Quantum physics with non-Hertian operators [^. It was found 
from numerical diagonalisation of the asymmetric BCS Hamiltonian that the transition 
between broken and unbroken PT-symmetry occurs on two boundary lines in the phase 
diagram. Furthermore, it was determined from the exact Bethe ansatz solution of the 
Hamiltonian that there exist free quasiparticle excitations on the boundary lines which 
can be characterised as satisfying generalised exclusion statistics. This finding relies 
on a deep understanding of the representation theory of a Yangian algebra, denoted 
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Y[gl{2)], from which the model may be obtained via the Quantum Inverse Scattering 
Method [?]. The objective of this article is to provide a detailed exposition of these 
generalised exclusion statistics. 

In Section 2 we begin by discussing the reduced BCS pairing Hamiltonian in general 
terms. We specify the asymmetric choice of Cooper pair scattering amplitudes which 
leads to our model of investigation. In Section 3 we review some relevant aspects of 
the Quantum Inverse Scattering Method |12lll9], and the associated techniques of the 
algebraic Bethe ansatz for the derivation of exact solutions. This algebraic approach is 
essentially an application of the representation theory for the Yangian algebra Y[gl{2)], 
which was developed by Chari and Pressley [6]. In Section 4 we illustrate how excitations 
obeying generalised exclusion statistics are identified in terms of dividing the solutions of 
the Bethe ansatz equations into two classes. This situation is reminiscent of generalised 
exclusuion statistics in an exactly solved model studied by Fendley and Schoutens [9], 
and provides an example of deconfined quantum criticality [TTIIIH] . In part our analysis 
amounts to addressing an old question regarding the completeness of the Bethe ansatz 
solution discussed in Bethe's original work for the XXX chain |5], and more recently 
debated in |2l[8]. We conjecture a formula for counting the number of deconfined states, 
and prove its validity in some limiting cases. Concluding remarks are given in Section 
5. Two Appendices are included. The first tables some complete, analytic solutions 
of the Bethe ansatz equations associated with low- dimensional Hilbert spaces of states. 
The second describes a mean-field analysis of the Hamiltonian which predicts that the 
Hamiltonian always has a real eigenspectrum. These latter calculations highlight the 
need for exact results to accurately characterise the system. 

2. The Hamiltonian 

The general form for a reduced BCS Hamiltonian as originally discussed in [3] is given 
by 

L L 

Hbcs = Yl ^^^j ~ Yl G'jfc4+4-s-s+- (1) 

j=i j,fc=i 

Here, j=l, . . . ,L labels a shell of doubly degenerate single particle energy levels with 
energies e^, and Uj = c^j^c^j^ + (^j-Cj_ is the fermion number operator for level j. The 
operators c-^, d^ are the annihilation and creation operators for fermions at level j. 
The labels ± refer to pairs of time-reversed states. 

An important feature of the Hamiltonian ([T]) is the blocking effect. For any unpaired 
fermion at level j the action of the pairing interaction is zero since only paired fermions 
are scattered. This means that basis states for the Hilbert space can be decoupled into 
products of paired and unpaired fermion states in which the action of the Hamiltonian 
on the space for the unpaired fermions is automatically diagonal in the basis of number 
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operator eigenstates. In view of this property the pair number operator 



L 



C"7'_L-C- ■ C"j 



commutes with ([T]) and thus provides a good quantum number. Below, M will be used 
to denote the eigenvalues of the pair number operator. 

It is convenient to express the Hamiltonian in terms of realisations of L copies of 
the su{2) algebra in the pseudo-spin representation, through the identification 

Si = c,_c,^, S^ = clc]_, S] = ^ (n, - J) . (2) 

The pseudo-spin operators satisfy the following su{2) commutation relations: 

[Sh S^] = ±6,,Sf, [SI, S,] = 26,,Si (3) 

Through this correspondence it is possible to identify a particular form of the 
Hamiltonian ([1]) as an exactly-solvable model associated with the su(2)-invariant six- 
vertex solution of the Yang-Baxter equation [7j. This occurs for the choice 



G 



jk 



G+ 

G+ + G- 


J <k 
j = k 
J >k 


2 



(4) 



for arbitrary G+ and G^. For real- valued G+ = G- this model corresponds to the case 
examined long ago by Richardson ^^U\. In [7] the specific case where G^ and G_ are 
a complex conjugate pair, leading to a self-adjoint Hamiltonian, was recognised as the 
Russian Doll BCS model of |T3] (see also [1]). 

Below we will instead consider the case for which G+ and G_ are both real- valued. 
Although this does not lead to a Hermitian Hamiltonian, nonetheless the spectrum is 
real- valued in some region of the coupling paramater space. Throughout we will work 
with a picket fence model whereby the e^ are uniformly and symmetrically distributed 
around zero. In particular we write 



J-^]S- (5) 



where the level spacing S provides an energy scale for the system. Through numerical 
diagonalisation of the Hamiltonian it is found that there are clearly identifiable regions 
where the spectrum is real. Illustrative cases are depicted in Fig. 1. It is evident that the 
boundary lines separating real-valued and complex-valued spectra are approximately 

G+-G^ = ±26 (6) 

for sufficiently large G+ and G_. 
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Figure 1: The shaded region depicts the values of the couphng parameters < G+, G_ < 
10 for which ( 1TM1I51) has real spectrum for the following choices: (a) 6 = 2, L = 8, m = 6, 
and M = 3; (b) (5 = 2, L = 8, m = 4, and M = 2; (c) (5 = 2, L = 8, m = 6, and M = 2, 
with blocked levels at 62 = —5 and eg = 3; {d) 6 = 2, L = 8, m = 6, and M = 2, with 
blocked levels at €4 = — 1 and €7 = 5. 



On the lines of the coupling parameter space given by ([6]) we will describe a class of 
free quasi-particle excitations exhibiting generalised exclusion statistics. Our analysis 
hereafter will be conducted using exact results. In order to present the exact Bethe 
ansatz solution for the Hamiltonian it is useful to make a change of variable. We 
parameterise the coupling constants G± through variables a, rj such that 



a 



2^nG: 



1 



V 



(G+ - G_) 



(7) 
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Inverting these relations gives 



_ 2r/e" 2r/e-" 

Gr_i_ — , Ct_ 



To simplify the discussion we will restrict attention to the subspace of unblocked states. 
The Hamiltonain is block-diagonal on this subspace with sectors determined by the 
eigenvalues M of the pair number operator N. For each M the dimension of the 
subspace within the space of unblocked states is L\/{M\{L — M)\). The total space 
of unblocked states has dimension 2^, with the exact solution for the energy spectrum 
given by 

M 

E = 2Y,v, (8) 

where the Vj are solutions of the Bethe ansatz equations [7] 

L M L M 

e'" n^^'^ - ^' + ^/2) n^^'^ - ^i - ^) = n^^'^ - ^^ - ^z^) n^^'^ - ^j- + ^)- ^9) 

1=1 j^k 1=1 j^k 

Before proceeding to an analysis of excitations characterised by the exact solution, 
we first need to recall and develop some results concerning the algebraic structure 
underlying the exact solution of the model. 

3. The 5f/(2)-Yangian, Quantum Inverse Scattering Method, and the 
algebraic Bethe ansatz 

The (7Z(2)-Yangian, which is denoted Y[gl{2)], is an infinite-dimensional algebra with 
generators {T^[l] '■ j,k = 1,2; I = 1, ..., oo}. In order to specify the algebraic relations, 
it is customary to use generating functions dependent on a complex variable [6j 

oo 

n{u) = Y,n[i]u-\ uec (10) 

1=0 

and then impose the commutation relations 

{u - v)[T,\u), Ttiv)] = v{Ttiv)Ti{u) - T^{u)Tl{v)). (11) 

For later use we note that the algebra admits a homomorphism : y[5f/(2)] — )■ ^[(7/(2)] 

ct>{Ti{u)) = p{u)Ti{u) (12) 

where p{u) is any function of u. 

Each irreducible, finite-dimensional, highest-weight y[5f/(2)]-module is charac- 
terised by a highest-weight vector |\E') such that [6] 

Tl{u)\m) = a{u)\^), 
Ti{u)\^) = d{u)\^), 
Tl{u)\^) = Q 
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where a{u) and d{u) are monic polynomials. Given such a module, we may obtain 
an equivalent module-action through f fT2l) with d{u) = p{u)a{u), d{u) = p{u)d{u). 
The condition on equivalences of module-actions is expressed through the Drinfeld 
polynomial P{u) [6] for which 



P{u — rj) a(u) 
P{u) ~ d(^' 

It is easily verified that 

a{u) d{u) 



(13) 



d{u) d{u) ' 
signifying the equivalence of the actions. 

Central elements are generated by the quantum determinant 

V{u) = Tl{u)Ti{u + r^)- Tl{u)Tl{u + r^) 
= Ti{u + r^)Tl{u) - Tl{u + r^)T',{u) 

which commutes with the generating functions fITOl) . The quantum determinant 'D{u) 
takes the eigenvalue n{u) = a{u)d{u + vi) on each irreducible, finite-dimensional, highest- 
weight module. In contrast to the Drinfeld polynomial the quantum determinant can 
distinguish equivalent module-actions, but might not distinguish inequivalent ones. 
For example consider a one-dimensional module with a{u) = d{u) = 1. Using the 
homomorphism flT2l) we then have for d{u) = d{u) = p{u) 

P{u) = 1, 

fi{u) = p{u)p{u + 1). 

Alternatively consider a highest- weight vector |\I') with the property that there exists 
f G C such that a{v) = d{v) = 0. It follows from the algebraic relations fITT]) that 

m = T^iv)m (14) 

satisfies 

Tlium = ~aiu)\^), 

with 

aiu) = a{uj, (15) 

u — V 

d(u) = "i^^^—ldiu), (16) 

u — V 

i.e. |\^) satisfies the conditions of a highest-weight vector. Here, either |\&) vanishes 
which is necessarily the case if the corresponding module is irreducible, or the module 
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with highest-weight vector |\E') is reducible. It is straightforward to check that |\1/) and 
|\^) admit the same eigenvalue for the quantum determinant T>{u), 

fi{u) = a{u)d{u + Tj) = a{u)d{u + tj). 

The fundamental two-dimensional y[(7/(2)]-module V{u) with basis {|1), |2)} 
admits the action 

Tl{u)\l) = («- V2)|l), Tl{u)\2) = {u + v/2)\2), 

T^{u)\l) = 0, T^{u)\2) = v\l), 

T^{u)\l) = v\2), T^iu)\2) = 0, 

T|(«)|l) = {u + r//2)|l), T|(«)|2) = {u - V2)|2). 

Since Y[gl{2)] is a bialgebra there is a co-product A : Y[gl{2)] -^ Y[gl{2)] ® Y[gl{2)] 
given by 

2 

A{Tliu)) = Y,Tn^)^niu-e) (17) 

1=1 

for arbitrary e G C. Iterating the co-product action permits the the construction of 
tensor product modules 

V{u; ei, ...cl) = V{u - d) ® .... ® V{u - e^) 

with highest weight vector 

1^) = |1)®^. (18) 

For generic choices of the ej this module is irreducible with Drinfeld polynomial [6] 

L 

P{u) = l[{u-e,+r^/2). (19) 

By using y [f7/(2)]-modules, abstract integrable models can be formulated as follows 
in the framework of the Quantum Inverse Scattering Method [T2l[T9]. For a G C define 
a transfer matrix as 

t{u) = e-°rii(M) + e°T|(M). 



From the commutation relations ( ITTjl it follows that 

[t{u), t{v)] = \fu,veC 

Through the above relation the transfer matrix is a generator of conserved (i.e. mutually 
commuting) operators for an abstract quantum system. A Hamiltonian can be defined 
as a polynomial function of the conserved operators. For the specific case of the 
Hamiltonian ([T]) subject to (jl]) we refer to [7]. 
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The algebraic Bethe ansatz provides a means to diagonalise the transfer matrix on a 
highest-weight module with highest weight vector |\1'). The eigenvectors of the transfer 
matrix are taken to be of the form 

M 

i=i 

where the ordering in the product is inconsequential, since the Tf{vj) commute. From 
the commutation relations flTTl) it is found that 



t{u)\Vi,V2,..;VM) = X{u)\Vi,V2,...,Vm) 

e~°'r] a{vk) ( jj Vk - Vj + r] 



with 



I I T-^{u)\Vi, ....,Vj-i,Vj+i...,VM) 

^ e"7/ d{vk) (^ Vk - Vj - v \ rT.2( ^ , \ 

Z.^73^|^n ,,_,^. JT,(n)K,.-,.,_„.,^,...,.M) 

(20) 



X{u) = e aiu) h e diu) 



U — Vj U — Vj 

J=l J -I 



Thus, if for all k = 1, ..., M, 



M M 

«(^fc) JJ {vk - Vj + r]) = e'^°'d{vk) JJ (vk - Vj - t]) , (21) 

the unwanted terms of the form \vi, ...., Wj_i, Vj+i...,vm) cancel. When this occurs either 
|fi,f2, ■.■,Vm) is an eigenstate of the transfer matrix, or |fi,f2, ■■■,Vm) vanishes. In the 
case that the state vanishes we refer to the corresponding solution set {t'i,f2, ...,fM} 
as spurious. In view of the homomorphism flT2l) . we can minimise the occurrence of 
spurious solutions sets in some instances by an appropriate choice of p{u) to express 
( 12T|) in the form 

M M 

P{vk -V)l[ K - Vj +v) = e^'^Pivk) n i^k - vj - 7]) . (22) 

For the choice of highest weight state given by ( IT8l) we obtain 

L 



\{{u - Q - r//2), 
1=1 

L 

d{u) = \[{u -ti + r//2) 



a[u) = 

1=1 

L 



1=1 
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and find that ( TT^ holds for (IT^ . Substituting the above expressions into (ETj) leads to 

Our main objective in the remainder of this work is to determine, for the specific 
case of the Hamiltonian ([1]) subject to (H|5|l and restricted to the hues IQ, those instances 
where |fi,W2, ■■■,vm) is an eigenstate and those instances where |fi,W2, ■■■,vm) vanishes. 
Through this investigation it will be seen how a class of free quasi-particle excitations 
exhibiting generalised exclusion statistics is uncovered. 

4. Generalised exclusion statistics 

Restricting to the lines (E]), which is equivalent to setting 6 = rj, the Bethe ansatz 
equations ([9]) are expressible as 

n(--('-i)) 

X (e'- (^v, + ^) f[{v, - V, -v)- (v, - ^) f[{v, - V, +v)\=0 (23) 

for each k = 1, ...,M. It is instructive to first consider the cases M = 1 and M = 2 
before addressing the general case. 

4.I. The one Cooper pair sector M = 1 

In this sector the dimension of the Hilbert space is L. Setting M = 1 and f i = f in fl23l) 
leads to 

L-l 



^^"l-f - "-f 



))n(^-.('-^))-. 



The above is a polynomial equation of order L, the roots of which can be stated explicitly. 
There is one root which is a-dependent, viz. 

r]L /e° + e-"~ 



^e" - e" 
while the remaining a-independent roots are elements from the set 

S = {6{k-L/2):k = l,...,L-l}. (24) 

To check whether any of these roots is spurious we consider the following action 
which is obtained through iterated use of the co-product (TT7|) : 

Ti^(«)i^)=^En(^-e.-v2) n(«-Q+v2)5;i^) 
j=i fc=i i=j+i 

L~l L 

= r^X[{u-5{k~L/2))Y,Stm. 
k=i j=i 
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For the a-dependent root, substitution into the above expression yield an eigenstate 
of the Hamiltonian. However, at first sight it appears that the roots from S are 
spurious. The state Tf(f)|\l/) vanishes for f G 5 due to the co-efficient polynomial 
11^=1 ("^ ~ '^(^ ~ -^/2))- Rescaling the state by this polynomial would then lead to the 
same eigenstate J2j=i^j~\'^) fo^ ^^^ v E S. However directly applying the Hamilonian 
to this state confirms that it is not an eigenstate. The problem with this approach is 
that rescaling by nA;=i (^ ^ ^(^ ~ -^Z^)) means the unwanted terms in (|20|) no longer 
cancel. 

Nonetheless we can confirm that the roots from S are not spurious. Starting with 
a generic value for rj we have 

L j-l L 

T^{u)\^) = r/ ^ J](n - e, - ^2) H {u - e, + v/2)S;\^) 
j=i k=i i=j+i 

fj.-l L 

For /i = 1, ..., L — 1 let 

fi-l L 

r,{u) = r]l[{u-e0-r]/2) J] (u-e^ + r]/2). 

/3=1 7=/^+l 

Then we define rescaled states 

U n'=,+i(^ - ^7 + V2) nfi(^ - e^ - V2) ^ ' ' ^ ' ' 

, y- m=l(^ - e. - V2) nt,+i(^ - Q + ^/2) ^ 

^ y- ng(» - 6^ - W2) nt,+i(^ - Q + ^/2) , 
,4tin^=,+i(^-e7+^/2)nS(«-^/^-W2) ^' ^' 

The Bethe ansatz equations ([9]) for M = 1 may be written as 

nU(^-e,-V2) ^ ^,, n;=i(^- ^7 + ^/2) 

n:?=^+i(« - ^7 + ^/2) nt,(^ - e. - r]/2) ■ 
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This then leads to 

{u-ep- r]/2) 



i^.(")>-En 'tr!!;f ^/i^>+^,^i*> 



2a y- U-e^ + ri/2 y^ U-ep + 7]/2 yr U - ej^i - T] /2 ^^ 

,4lti ^' - ^^- - ^/2 ii « - ^/^+i + ^/2 i,ii ^ - ^^ - ^/2 ' 

Next let ?7 = (5 giving 

|$,(5(;. - L/2))) = ^ 5, |vl/) + e ^>_^^ (^-L/2)^-(L- 1)^2-^2 ^^- '^^ 

1^ 2a ^ 

=E^;i^)+^E^/i^)- (25) 

It can be verified by direct calculation that the expressions for |$^(5(yU — L/2))), /x = 



1, ...,L — 1 as given by ( 1251) are eigenstates of the Hamiltonian ([T]) subject to (I4ll5]l6l) . 
Along with the eigenstate associated with the a-dependent root, this provides a complete 
set of L eigenstates for the sector M = 1. 

4-2. The two Cooper pair sector M = 2 

In this sector the dimension of the Hilbert space is L{L — l)/2. The first important 
observation to make from the previous Subsection is that when r] = 6 the roots of 
the Bethe equations can be clearly demarcated into two groups, those which are a- 
independent and those which are not. This provides three sub-cases to consider. 

4-2.1. Two a -independent roots. From the Bethe ansatz equations (!23|) it is seen that 
a solution set is formally obtained by choosing any two elements from the set S given 
by Eq. fl24|) . including the case when the roots are equal. In this respect the excitations 
have the character of two free quasi-particles. However some solution sets are spurious. 
To identify those instances, we first note that the one Cooper pair sector eigenstates (!25|) 
are F [(7/(2)] highest weight states of the form (fT4l) with suitable rescaling. Specifically, 
choosing vi = 6{fi — L/2) leads us to consider 

Tl{u)\^,i6ifi - L/2))) = ~aiu)\^Ml^ - L/2))), 
r|(«)|$^(5(/i - L/2))) = d{u)\<l>Mf^ - L/2))) 
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with 

L 

a{u) = {u-S{fi-l- L/2)) Y[{u - S{1 - L/2)), 

L 

d{u) = {u-5{ii + l- L/2)) Yl {u-6{l-l- L/2)) 
and the associated Drinfeld polynomial 

L 

P{u) = Y[ {u-6{l-l- L/2)) 

such that flT3l) holds. For this Drinfeld polynomial we find that the set of a-independent 
roots for the Bethe ansatz equations (I22|) . from which we can choose f2, is the restricted 
set 

5' = {5{k~L/2) : k = 1,...,L-1; fc 7^ /x - l,/i,/x + 1}. 

Here we observe the manifestation of generalised exclusion statistics. Having first chosen 
the Bethe root vi = 6{fi — L/2) we find that not only can we not choose it again (as 
is the case for the familiar fermionic exclusion principle), we also cannot choose the 
"neighbouring" roots 5(/i± 1 — L/2). Bearing this in mind, a simple counting argument 
shows that the number of eigenstates in the M = 2 sector where both roots are a- 
independent is given by (L — 2)(L — 3)/2. 

4.2.2. One a-independent root. Choosing vi E S the Bethe ansatz equation for V2 ^ S 
is the quadratic equation 

e2"(t;2 + 6L/2){v2 ~ v^ - S) = iv2 - SL/2){v2 - vi + S). 

The solution reads 

_ 2(1 - e^")t;i + ^(1 + e^")(L -2)±Vd 
^'~ 4(1 - e2°) 

where 

D = 6^6^'' + 1)\L - 2)2 + 4(e2" - l)\vf + 26^L) + 45(e^" -1){L + 2)vi. 

While there are generally two solutions for V2, two special cases need to be re- 
examined in closer detail. Choosing vi = (5(1 — L/2) the equation for V2 is 

e2"(t;2 + 5L/2){v2 - v^ - 5) = iv2 - 5L/2){v2 - Vi + 5) 
e^''{v2 + ^L/2){v2 - 25 + bL/2) = {v2 - 6L/2){v2 + 5L/2). 
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Formally there are two solutions for V2, but the case V2 = —SL/2 is a- independent and 
spurious by the reasoning presented in the preceeding Subsubsection. The a-dependent 

root is 

(L/2 + (L/2 - 2)e2")(5 
-2 = Y^^. • 

Similarly, choosing vi = 6{L/2 — 1) the equation for V2 is 

e2°(t;2 + 5L/2)iv2 -v^-5) = (t;2 - 6L/2)iv2 - v^ + 5) 
e'^'^{v2 + 6L/2){v2 - 6L/2) = {v2 - 6L/2){v2 - 6L/2 + 26). 

In this instance the a-independent solution t>2 = 6L/2 is spurious. The a-dependent 

root is 

((L/2 - 2) + e2"L/2)5 

'' = T^^^ • 

These results indicate that the number of eigenstates in the M = 2 sector, where only 
one root is a- independent, is given by 2 (L — 2). 

4-2.3. No a-independent roots. Since the the dimension of the Hilbert space for M = 2 
is L{L — l)/2, and we have accounted for (L— 2)(L— 3)/2 solutions for two a-independent 
roots and 2(L — 2) solutions for one a-independent roots, there can only be one non- 
spurious solution with both roots being a-dependent. Although we have not been able to 
derive this result in the sense of producing an explicit general formula for the solution, 
we have checked the low-dimensional cases L = 4, 5, 6 with the results tabulated in 
Appendix A. Specifically, we have verified that for these cases our analyses above are in 
complete agreement with results obtained by direct diagonalisation of the Hamiltonian 
to obtain the eigenspectrum. In particular, direct diagonalisation confirms the picture 
that the a-independent roots are associated with free quasi-particle excitations with 
generalised exclusion statistics. 

4.3. The case of general M 

For a given L and M let P, with P < M, denote the number of a-independent 
roots within a set of M roots. Further let n{L, M, P) denote the number of non- 
spurious solutions for each set of these quantities. From our previous discussions we 
have n{L, 0, 0) = 1, n{L, 1, 1) = L - 1, n(L, 1, 0) = 1, n{L, 2, 2) = (L - 2)(L - 3)/2, 
n{L, 2, 1) = 2(L - 2) and n(L, 2, 0) = 1. 

Next we will prove a formula for n{L,M,M), i.e., the case when all roots are 
a-independent. We saw in the previous Subsection for M = 2 that if we choose a 
particular root Vi = 6{fi — L/2), the generalised exclusion principle prohibits the choice 
f 2 = 6{fi ± 1 — L/2) to obtain an eigenstate. This result generalises for arbitrary M in 
a straightforward manner, viz. given a solution set for the M — 1 Cooper pair sector 

{vj = S{fi,-L/2):j = l,...,M-l} 
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we are prohibited from choosing any neighbouring roots S{fij ± 1 — L/2) for vm to 
obtain an eigenstate in the M Cooper pair sector. The proof is a matter of iterating the 
procedure described for M = 2 in Subsubsection 14.2. II To determine n{L, M, M) we 
need to count the number of a-independent solution sets which respect to this exclusion 
rule. 
Let 

P 

Q 

denotes the number of ways that q identical objects can be placed in p boxes such that 
the objects cannot be placed into adjacent boxes. This quantity satisfies the recursion 
relation 

p — 1 



p 



+ 



p-2 
q-1 



with the initial conditions 



p 




q 

1 
1 



p>0, 
q>0, 



1. 



The recursion relation can be solved to obtain 



p 



(p + l-g)! 



{p + l-2q)\q\ 
The solution set S as given by f l24|) contains L — 1 elements, yielding 



n{L,M,M) 



L-1 

M 



{L-M)\ 
{L-2M)\M\' 



Although we have not been able to prove a general result for n[L,M,P), we 
conjecture that it is given by 

{L-M)\M\ 



n{L,M,P) 



(26) 



(L-M-P)!(Af-P)!(P!)2' 

This formula agrees with the subcases P = M, and M = 0, 1, 2 with P < M, discussed 
above. It also satisfies the required condition 

M 

J2n{L,M,P) 



L\ 



p=o 



{L-M)\M\ 



to account for all states in each sector of M Cooper pairs. The above identity may be 
proved by considering binomial expansions of the terms in 
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5. Conclusion 

The Hamiltonian ([T]), subject to ( Hll5l) . admits an exact solution. On the subspace of 
unblocked states the energy eigenvalues are given by ([8]) where the set {vi, ...vm} is a 
solution of the Bethe ansatz equations ([9]). The energy expression is a simple sum of the 
Bethe roots, which leads to a natural interpretation of each Bethe root being associated 
with a quasi-particle. In this quasi-particle picture the states are generally bound rather 
than free, since a solution set for M quasiparticles is not simply the union of one-body 
solutions due to the coupled nature of ([9]). In general the roots Vj will depend on the 
two coupling parameters G+, G_, or equivalently, the variables a and r] as given by (^^. 

Setting rj = 6 corresponding to ([6]), where S is the level spacing in (|5]), we find that 
the roots of the Bethe ansatz equations divide into two classes, viz, those which are a- 
dependent and those which are a-independent. From this perspective we made identify 
the lines (Q with deconfined excitations, those being associated to the a-independent 
roots which do not occur for general 77. From the form of the Bethe ansatz equations 
given by fl23!l . the a-independent roots were determined explicilty. They are elements 
of the set S defined in ( !24l) . The corresponding states were found to be highest weight 
states of the underlying Y[gl{2)] algebraic structure. These roots are associated with 
free quasi-particles excitations, since a solution set for M quasiparticles is a union of 
one-body solutions. However, a close examination of the states associated with these 
solutions led to the conclusion that there are spurious solution sets, which give rise to 
an interpretation of generalised exclusion statistics. We have conjectured the formula 
fl26|) for the number of non-spurious solution sets with M roots, where P of the roots 
are a-independent. A proof of the result for P = M was provided, and it was also found 
that ( l26l) is valid in certain limiting cases. 

It was proposed in [TTlfTS] that quantum criticality may be identified by deconfined 
excitations at the critical point, which are not found in phases adjacent to the critical 
point. It is interesting to note for the present study that the lines ([6]) on which 
the deconfined a-independent excitations occur are in very close agreement with the 
boundary lines between the unbroken PT-symmetric phase and a broken "PT-symmetric 
phase. The PT-symmetry breaking phase boundary lines are clearly identified in Fig. 
1, which was obtained by direct numerical diagonalisation of the Hamiltonian. 

Appendix A - Tables of roots of the Bethe ansatz equations 

In this Appendix we list the roots of the Bethe ansatz equations for some small sized 
systems L = 4, 5, 6, M = 2 and no blocked states, for when the coupling parameters 
correspond to the lines (I6l). The first column lists the energy eigenvalues obtained 
by direct diagonalisation, while the second column lists the roots of the Bethe ansatz 
equations ( 123|) associated to each energy eigenvalue through (IHl). These results are 
consistent with conjecture (!26|) concerning the distribution of roots, and establish that 
the energy spectrum of the Hamiltonian on the lines ([6]) is real-valued for these cases. 
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E 


(^1,^2) 





(5, -5) 






(1 + e2° 


l_e2a 


/ (1 + e"" + V9 - 14e2" + 9e'^^)6\ 
V ' 2(1 - e2") ) 


(1 + e^" 


- V9 - 14e2" + 9e4°)5 
l-e2" 


/ (1 + e^" - V9 - 14e2° + 9e4°)5\ 
V ' 2(1 - e2") J 


l_e2a 


/ (3 + 3e2" + VI - 14e2" + e^") 5 
1^ 2(l-e2") 


(3 + 362^^ - VI - 14e2" + e^") 5 \ 
2 (1 - e2°) J 



Table 1: The energy spectrum and associated solutions of the Bethe ansatz equations 
fl25]) on the lines (^ ioi L = A, M = 2. While the pairs (5,0), (0,-5) are solutions 
of (123!) they do not correspond to an eigenvalue of the Hamiltonian, giving rise to a 
generalised exclusion principle. 



Appendix B - Mean field analysis 

In the Appendix we show that a mean-field annalyis in the continuum limit gives the 
result that the spectrum of the Hamiltonian ([1]) subject to ( Hll5l) is real for all values of 
the coupling parameters. For this reason it is deemed necessary to analyse the model 
through the exact Bethe ansatz solution. 

For a general non-Hermitian Hamiltonian with a real spectrum, where |\1/^) denotes 
the right ground-state eigenvector and (^'| denotes the left ground-state eigenvector, the 
ground-state energy can be expressed as E = (^^|i7|^'')/(\&^|\E''') provided (\1/'|\E''') 7^ 0. 
When this is the case, for any operator A we define the ground-state expectation value 
as (A) = (\E'^|74|\l/'')/(\l/^|\l'''). This definition preserves the Hellmann-Feynman theorem. 

Given a BCS Hamiltonian of the form (Q, we introduce a (real- valued) chemical 
potential /x and set 



Aj- = ^Gjk{ck-Ck+) 



^l 



E^. 



jk 






k=l 
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E 



1^1, ^^2j 



26 



36 6 

T'~2 



3(5 _36 



-26 



'6_ _36^ 
2' 2 



2 (2 + e^") 5 



'35 (l + 5e2"),5^ 



1-e 



2a 



2 ' 2(1 -e' 



2q^ 



2 (1 + 2e2") ,5 



1-e 



2a 



35 (5 + e^") 5 
"T'2(l -e 



2a ^ 



(3 - V9 - 16e2« + 16e4") 5 



1-e 



2a 



5 (2 + e^" - V9 - 16e2° + IGe^") 6 
2' 



2 1-e 



2a ^ 



(3 + V9 - 16e2« + 16e4°) 5 



1-e 



2a 



6 (2 + e2° + V9 - 16e2" + IGe^") 6 
2' 



2 1 - e 



2a ^ 



(3e2" - V16 - 16e2° + Oe^") 5 



1-e 



2a 



5 (l + 2e2" - V16 - 16e2" + Qe^") 5 
"2' 



2 (1 - e' 



2a A 



(3e2" + V16 - 16e2" + Oe^") 5 



1-e 



2a 



5 (1 + 2e2" + Vl6 - 16e2" + <de^^)6 

'2' 



2 1-e' 



2a ^ 



(4 + 4e2" + Vl - 18e2° + e^") 5 



:i + e^") 5 



2(1 -e 



2a\ 



1-e 



2a 



(4 + 4e2" - Vl - 18e2" + e^^") 5 



2(1 -e 



2a ^ 



Table 2: The energy spectrum and associated solutions of the Bethe ansatz equations 
(I23D on the lines ([6]) for L = 5, M = 2. While the pairs (35/2,5/2), (5/2,-5/2), 
(—5/2, —35/2) are solutions of f l23|) they do not correspond to an eigenvalue of the 
Hamiltonian, giving rise to a generalised exclusion principle. 



to obtain the mean-field approximation 

L L L L 

HuF = Y^ ejHj - Y^ A^.c]+c]_ - Y Kck-Ck+ + ^ Gjk (b]j (6^) -fi{n-{n)) 

j=l j=l fc=l j,k=l 



Setting ^fc = ek — fi, £k = a/^| + A^A^, by diagonalising i^MF it is seen that the 
elementary excitation energies are simply the Sk- The entire energy spectrum is real- 
valued if the ground-state energy is real-valued and the products Aj^ A^ are non-negative 
for all k. 
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E 



(^^1,^2) 



46 



(25,0) 



26 



(26,-6) 



(6,-6), (26,-26) 



-26 



(6, -26) 



-46 



(0, -26) 



(5 - e^" - V9 - 18e2" + 25e4") 6 



(3 + e^" - V9 - 18e2° + 25e4") 5 



1-e 



2a 



2('l-e 



2a 1 



(5 - e^" + V9 - 18e2" + 25e4") 5 



1-e' 



2a 



(3 + e^" + V9 - 18e2" + 25e4") 5 
2(l-e2") 



2 (1 + e2" + 2V1 - e2" + e^") 5 



1-e' 



2a 



0, 



(1 + e2" + 2V1 - e2" + e4") 5 



1-e' 



2a 



2 (1 + e2" - 2V1 - e2° + e^") 5 



(1 + e^" - 2V1 - e2« + e4°) 5 



1-e' 



2a 



1-e' 



2a 



(5e2" - 1 - V25 - 18e2" + 9e*») 6 



(1 + 3e2° - V25 - 18e2" + 9e4") 5 



1-e' 



2a 



2 (1 - e' 



2a ^ 



(5e2° - 1 + V25 - 18e2" + Qe^") 5 



1-e' 



2a 



(1 + 3e2° + V25 - 18e2« + 9e4")(5 



2(1 -e' 



2a\ 



2 (3 + e^") 5 



1-e' 



2a 



26, 



[1 + 3e2")(5 



1-e 



2a 



2(l + 3e2")5 



1-e' 



2a 



-26 



(3 + e^°) 6 
' 1 - e2" 



10 (1 + e^") 5 



1-e- 



2a 



(5 + 5e2" + VI - 22e2" + e^") 5 
2(1 - e2") ' 

(5 + 5e2" - VI - 22e2° + e^") 5' 



2 (1-e 



2a\ 



Table 3: The energy spectrum and associated solutions of the Bethe ansatz equations 
(I23D on the lines (E]) for L = 6, M = 2. While the pairs (26, 6), (6, 0), (0, -6), (-6, -26) 
are solutions of fl23l) they do not correspond to an eigenvalue of the Hamiltonian, giving 
rise to a generalised exclusion principle. 
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The right mean-field ground state is given by 

L 
k=l 

where vl/ul = {Sk — ^k)/^l = ^kli^k + ik)- Analogously for the left ground state, 

L 

(vp'i=(oin(«u+4cfc-cfc+), 



fc=l 



we have v\/u\ = {£k — ^k)/'^\ = ^\/{£k + ik)- Then (\l/i|\l/i.) is found to be non-zero 
provided A),A^ ^ for all k. 

Self-consistency requirements impose that 



fc=l ^ 

L 

) = L + J2^- (29) 



^l = lj2^^'^% (28) 



We refer to fl27|28p as the as the gap equations, and to fl29l) as the chemical potential 
equation. Using these equations allows for the mean-field ground-state energy to be 
expressed as 



fe=l ^ '^^ fc=l '^ 



Note that for a Hermitian Hamiltonian of the form ([T]), A\, A^ are a complex 
conjugate pair, in which case ( I27f28!) are equivalent. The question remains whether 
for non-Hermitian Hamiltonians there exist solutions of (127112811291) such that AlAI is 



^k^k 

real-valued for all k. For the choice (jl]) we next show that this is the case in the 
continuum limit. We introduce a cut-off energy u by setting 6 = 2u/{L — 1) such that 
ei = —u and e^ = u. Letting G±. = 4ujg±/L, x = (n)/(2L), in the continuum limit 
L — )• oo, G± — !■ 0, (5 — )■ de' we have fl27f28ll29l) assuming the integral equation forms 
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where £{e) = y/{e - iJ,y + A\e)A'^{e). Differentiating (^ and (EH) leads to the 
conclusion that A = ^yA^(e)A^(e) is constant. With this observation the integrals 
(130)13111321) can be evaluated to obtain 



_ ujgl + g'^){2x - I) ^ _ lQ(^^glg-x{l - x) 



gl-g- ' {gl - g-)^ 

where x = V(fl'+ ~ 5'-)- Note that A^ > for all x, g±. Now the elementary excitation 
spectrum is explicit, viz. £{e) = iy(7^^70^~+"^! ~^ < e < i^, such that A is the gap. 
We obtain the gound-state energy per fermion as 

,. Emf 
cmf = hm — — 

L^oo ZxL 

1 /•" ^^ 2e(e -fi) + A^ 



6xu 
1 



ixoj 



v/(e - fxy + A2 

(w + fi)^/iuj - fiy + A2 + (w - i2)y/{uj + fiy + A2) . 



Thus within this mean-field analysis the non-Hermitian Hamiltonian has a real spectrum 
for all couplings g± and filling fractions x. In the Hermitian limit g± -^ g for which 
X -^ C)0, both fi and A can be evaluated through use of exp(a;) = lim (1 + x/n)"'. In 
particular for half-filling x = 1/2 we obtain /x = and A = u;/sinh(l/25f) which is in 
agreement with the classic s-wave result obtained in [3] (equation (2.40)). 

It might be expected that mean-field results are exact in the thermodynamic limit 
(e.g. see [16] for when g+ = gJ). It would be very useful to determine whether or not the 
energy spectrum is real to leading order, with complex terms only appearing in lower 
order corrections. 
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